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BORIS VERTMAN 

Abstract. We consider a combinatorial Laplace operator on a sequence of discrete 
graphs which approximates the m-dimensional torus when the discretization pa¬ 
rameter tends to infinity. We establish a polyhomogeneous expansion of the re¬ 
solvent trace for the family of discrete graphs, jointly in the resolvent and the dis¬ 
cretization parameter. Based on a result about interchanging regularized limits and 
regularized integrals, we compare the regularized limit of the log-determinants of 
the combinatorial Laplacian on the sequence of discrete graphs with the logarithm 
of the zeta determinant for the Laplace Beltrami operator on the m-dimensional 
torus. In a similar manner we may apply our method to compare the product of the 
first N C N non-zero eigenvalues of the Laplacian on a torus (or any other smooth 
manifold with an explicitly known spectrum) with the zeta-regularized determinant 
of the Laplacian in the regularized limit as N —> oo. 
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1. Introduction and formulation of the results 

Introduction of a zeta-regularized determinant for the Hodge Laplacian on com¬ 
pact Riemannian manifolds by Ray and Singer [RAS171] provided a counterpart to 
the determinant of discrete Laplacian on a simplicial complex associated to a tri¬ 
angulation of the manifold. Relation between the zeta-regularized determinant of 
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the Hodge Laplacian as an analytic spectral invariant, and the determinant of the 
discrete Laplacian as a combinatorial quantity, has been shown to go beyond be¬ 
ing just formal counterparts by the proof of the Ray-Singer conjecture [RAS171] by 
Cheeger [CHE79] and Muller [MUL78]. 

In fact, Muller [MUL78] proved that a specific combination of determinants for 
discrete Laplacians in various degrees, defined on simplicial complexes associated 
to a triangulation of a compact Riemannian manifold, converges to the correspond¬ 
ing combination of zeta-regularized determinants of Hodge Laplacians when the 
mesh of the triangulation goes to zero. Another instance of a link between determi¬ 
nants of discrete Laplacians and zeta-regularized determinants of the correspond¬ 
ing Hodge Laplacians is a recent joint work with Reshetikhin [REVE13] which in part 
motivated the Burghelea Friedlander Kappeler gluing formula for determinants by 
studying the corresponding combinatorial problem. 

In both instances the behavior of the individual determinants remained open, 
since the discussion is rather based on existence of a well-defined limit for combina¬ 
tions of determinants for discrete Laplacians under finer discretizations. This leads 
to the general problem if the zeta-regularized determinant of a Hodge Laplacian 
may indeed be recovered from its discretization. This translates into a question on 
existence of an asymptotic expansion for determinants of discrete Laplacians under 
refinement of discretization. 

Interest in the asymptotic behaviour for determinants of discrete Laplacians arises 
in various mathematical settings, even without the conjectured relationship with 
its zeta-regularized counterpart. In statistical mechanics the interest stems from 
identification of the determinant in terms of the number of spanning trees on a 
graph by Kirchhoff [K1R47]. Moreover, in the setting of two-dimensional lattices, 
determinants of certain Z 2 subgraphs were expressed in terms of the number of 
dimer coverings of related Z 2 subgraphs by Temperley [TEM74]. 

In mathematical physics, existence of an asymptotic expansion for determinants 
of discrete Laplacians may provide a way to construct quantum field theory of a free 
scalar Bose field as a scaling limit of a Gaussian quantum field theory on a discrete 
simplicial complex associated to a triangulation of the manifold, as the mesh of the 
triangulation goes to zero. In fact, this intuition also lies behind Hawking [HAW77]. 

In fact, asymptotics for determinants of discrete Laplacians has been studied in 
several instances with partial results. In the setup of rectilinear polygonal domains, 
Kenyon [Kenoo] derived a partial asymptotic expansion for the determinant of the 
corresponding discrete Laplacian. However existence of a constant term in that par¬ 
tial expansion, let alone its identification with the zeta-regularized determinant, re¬ 
mained an open problem. Other related results include Burton-Pemantle [BUPE93] 
and Sridhar [SRI13]. 

In the setting of tori, the spectrum of the discrete and Hodge Laplacians is under¬ 
stood explicitly, which naturally allows for finer asymptotic results. This setting has 
been studied by Chinta, Jorgenson and Karlsson [CJK10], who equated the constant 
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term in the asymptotics for determinant of discrete Laplacians with the logarithm 
of the zeta-determinant, cf. also the preceeding results by Kasteleyn [Kas6i] and 
Duplantier-David [DuDa88]. Their analysis is based on a discussion of the discrete 
heat operator in terms of Bessel functions, and is strongly rooted in the explicit 
structure of the spectrum. In the setting of two-dimensional tori, Chamaud [Chao6] 
has relaxed geometric assumptions by variational methods. 


Closely related to the question, if zeta-regularized determinant of a Hodge Lapla- 
cian may be recovered from asymptotics of its discrete counterpart, is a problem of 
relating the zeta-regularized determinant to finite eigenvalue products. Asymp¬ 
totic behavior of eigenvalue products has been studied by Szego [SZE15] for certain 
Topelitz matrices, and in fact Friedlander-Guillemin [FrGuo8] have compared the 
Szego and zeta-regularized determinants for zero-th order pseudo-differential op¬ 
erators. In a related work Friedlander [FR189] obtained the zeta-determinant of a 
higher order elliptic pseudo-differential operator A by considering the asymptotics 
of a determinant for some determinant class operator associated to A 1 . 


The presented references share the common idea of replacing the meromorphic 
continuation technique used in the definition of the zeta-determinant, by analysis 
of asymptotic expansions of classical determinants. Our paper studies this ques¬ 
tion in the explicit setting of m-dimensional tori and reproves the result Chinta, 
Jorgenson and Karlsson proposing an alternative independent ansatz. Our method 
is based on a polyhomogeneous asymptotic expansion of the combinatorial resol¬ 
vent trace jointly in ther resolvent and the discretization parameters. In particular 
this approach may be be viewed as part of a program initiated jointly with Lesch 
[LEVE13], cf. also Vertman [VER13] and Sauer [SAU13]. The main technical tool is 
a careful analysis of the terms in the Euler Maclaurin formula, as well as a result 
on interchangeability of regularized limits and integrals. We focus on identification 
of the regularized limit of discrete determinants in terms of the zeta-regularized 
determinant, rather than studying other terms in the asymptotics as in [CJK10]. 


Remark 1.1. Our discussion is applicable beyond the setting of discrete tori. A brief 
look into our argument makes apparent that our results depend on existence of 
a polyhomogeneous expansion of the combinatorial resolvent trace, which needs 
not be an exclusive feature of m-dimensional tori. We therefore expect our results 
to have applications to explicitly computable quotients of M m under action of Z m 
lattice subgroups. 

Remark 1.2. We point out that the same principle may be applied to identify the 
product of the first N e N non-zero eigenvalues of the Laplacian on a torus (or any 
other smooth manifold with an explicitly known spectrum) with its zeta-regularized 
determinant in the regularized limit as N —> 00. 
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1.1. Hadamard partie finie regularization. Consider f e C°°(M+,C),M + := ( 0 , oo) 

such that for x —4 oo in the Landau notation 

N-1 Mj M 0 

(l.l) f(x) = Y_ Y a fr X<Xi log^M + Y a 0klog k (x) + o(x aN log MN (x)) 

j=l k=0 k=0 

for some N e N and (<Xj) C C, such that (Re(cq)) is a monotonously decreasing 
sequence with Re(ctN) < 0 . Then we define the regularized limit of f(x) as x —> oo by 

LIMf(x) := a 0 o. 

x —>00 

If Re(aN) < — 1 , the integral of f over [ 1 , R] admits an asymptotic expansion of the 
form (1.1) as R —» oo, and we set 


- f(x)dx:=LIM 

R—>oo 


f(x)dx. 


Similar definition holds for the regularized limit at x = 0 assuming an appro¬ 
priate asymptotic expansion of f(x) as x —> 0 , of the form (1.1), where (Re(cq)) 
monotonously increasing with Re(<XN) > 0. 

One crucial analytic property of regularized limits and integrals is the following 
interchangeability result, which has been presented in [LeVei3, Lemma 3.3] under 
mildly stronger assumptions on the asymptotics 1 . We provide a proof for general 
asymptotic expansions of the form (1.1). 

Proposition 1.3. Let f e C°°(R^,C),M + = ( 0 , 00), be homogeneous of order d e C 
jointly in both variables and asymptotic expansions of the form (1.1) in each of the variables 
(z, n) —> 00 individually (with the other variable fixed) 

N-1 Mj M 0 

f(z» 1) = Y Y a ’ kZtXi lo g k ( Z ) + Y a 0 klog k (z) + o(z“ N log MN (z)), 

j=1 k=0 k=0 


N1 


M' 


f(l,n) = Y Y b ) knPi lo g k ( n ) + Y- ^oklog k (n) + o(n PN ' log M N(n)), 

j=1 k=0 k=0 

for some N,N'gN and (cq), (|3j) C C, such that (Re(cq)), (Re((3j)) are monotonously 
decreasing sequences with Re(a N ) < — l,Re(( 3 N /) < min{ 0 , d+ 1 }. Then 


(1.2) 


LIM - f(z, n)dz= - LIMf(z, n)dz + Corr, 


where Corr = f (z, 1 )dz if d = — 1 and zero otherwise. 

It is an inherent part of the statement, that the regularized limits and integrals in 
the equality (2.2) exist. 


1 In contrast to [LEVE13, Lemma 3.3] we do not assume that f(z, 1 ) and f( 1 ,n) are smooth at 
z, n = 0 , and require partial asymptotics as z, n —1 00 instead. In fact the latter result is applied to 
some homogeneous f(z, n), where smoothness of f(z, 1) and f (1, n) at z, n = 0 indeed fails. 
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1.2. Polyhomogeneous expansion of combinatorial resolvent traces. For any in¬ 
teger n G N we consider the quotient space Z/uZ which we refer to as a discrete 
circle. A finite product of m copies of Z/nZ defines a discrete torus T™. which may 
be viewed as a discretization of the m-dimensional torus manifold T m , given by a 
product of m copies of S 1 . Here, S 1 C l 2 is a circle with radius 1 . 


The combinatorial Laplacian A n on the discrete torus T™ is the sum of the Lapla- 
cians L n on each discrete circle Z/nZ component, defined for any f : Z/nZ —> Mby 
the following difference operator 

OCnfHM) = ^2 (( f (M) - f([k -1])) + (f(M) - f(Dc +1]))) • 

Then the spectra of £ n and A n amount to 


(i- 3 ) 





k G N 0 , k < u 


} 



lq G N 0 , lq < n for i = 1 ,..., 



In this multiset notation, eigenvalues A appear multiple times according to their 
multiplicity m(A). By the particular choice of the rescaling factor ^ 4 , eigenvalues 
of L n and A n approximate the eigenvalues of the Laplace Beltrami operators on S 1 
and T m respectively, as n —> oo. 


Consider the combinatorial Laplacian A n on the discrete torus T™ and traces of 
the corresponding resolvent powers Tr(A n + z 2 )“ for any a G N and z G M + . We 
refer to the quantities Tr(A n +z 2 ) _m and TrfA+z 2 ) m as resolvent traces. The resolvent 
trace of the Laplace Beltrami operator A on the smooth torus T ,Tl can be obtained as 
the limit of the combinatorial trace. In fact we have the following result. 


Proposition 1.4. The resolvent traces are smooth in z G M + and for k G N 


lim 3 ' k) Tr(A n + z 2 


0 ' k) Tr(A + z 2 r m . 


Our first central result is a polyhomogeneous expansion of the resolvent trace for 
the combinatorial Laplacian A n , jointly in the resolvent parameter z G M+ and the 
discretization parameter n G N. 

Theorem 1.5. The resolvent trace for A n admits a partial polyhomogeneous expansion 

m 

Tr(A n + z 2 )“ m = Y h_ m -j(z,n) + H(z,n), 

i=o 

where each h_ m _j G C°°(M 2 ) is homogeneous of order (—m—j) jointly in (z, n), h_ m _j (z, 1 ) 
and Tx._m._j C1 ,admit an asymptotic expansion of the form (1.1) as z,n-) 00, respectively. 
The remainder term satisfies H N (z, n) = 0 (z~ 2m ~ 2 ), as z —> 00, uniformly in n > 0 . 
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1.3. Zeta-functions and zeta-regularized determinants. In the next step we intro¬ 
duce (zeta-regularized) determinants of A n and A. In the discrete case, the deter¬ 
minant of A n is defined here as a product of its non-zero eigenvalues, counted with 
their multiplicities, and in fact satisfies the following integral representation 


(i- 4 ) 


log det A n 


-2 


-- zTr(A n + z 2 ) 'dz, 


Jo 

which is an immediate consequence of the following computation 


^°° zdz 

o (A + £ 2 ) 


-LIMLIM 

R —^00 g—>0 


log (A + 



z=e 


log(A), if A ^ 0 , 
0 , if A = 0 . 


Integrating (1.4) by parts (m — 1 ) times (more precisely we perform integration by 
parts for the integral on [e, R] and take the regularized limit of the total expression 
as e —> 0 and R —> 00) yields 


(i- 5 ) 


log det A n = —2 j- z 2m l Tr (A 


z 2 r m 


dz, 


with the boundary terms vanishing in the regularized limit. 


The zeta-regularized determinant of the Laplace Beltrami operator A is obtained 
by the following procedure. The zeta-function of A is defined for Re(s) > m /2 by 


C(s, A) := Y_ m(A)A s , 

AeSpecA\{0} 


where m(A) denotes the multiplicity of the eigenvalue A. The integral expression 
(cf. [LEVE13, Section 1.3]) amounts after iterative integration by parts to 


C(S, A) =2 


sin7ts P(1 — s)r(m) 
7T r(m— s) 


-- Z' 


2m—2s—1 


Tr(A + z 2 )^ m dz, 


with the standard asymptotic expansion of the resolvent trace Trf A + z 2 ) 111 yields a 
meromorphic extension of a s , A) to the whole complex plane C with s = 0 being a 
regular point. We define the zeta-regularized determinant by 


(1.6) 


log deh A := — 


ds 


s=0 


as, A) = -2 f z 2m_1 Tr(A + z 2 r m dz 


,2',—m. 


1.4. Main result: Approximation by combinatorial determinants. 

Our main result now reads as follows. 

Theorem 1.6. The logarithmic determinant log det A n admits a regularized limit as n —> 
00, which equals the logarithm of the zeta-determinant of the Laplace Beltrami operator, i.e. 

log dek A = LIM log det A n . 

n—>00 
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A remark on the relation to the discussion in [LEVE13] is in order. By Proposition 
1.4 we may write 


log det c A = -2 4 - z 2m - 1 Tr(A n + z 2 r m dz 




0 

poo 


2m—1 

0 n— >00 

■*00 


= -2 f z Zm ~' lim Tr(A n + zT m dz 

n—1 n—1 


= —2 — z 2m 1 lim y ■ ■ ■ y (co(n,k!, ...,k m ) + z 2 ) m dz, 

A TL—>00 LL 7 

•'U V_V —C\ 


ki —0 k m —0 


where we have introduced 


(i- 7 ) 


cu(n, ki,k m ) := ^sin 2 ( ^ j , 


i=1 


Tlk; 


and hence our main Theorem 1.6 looks as an application of the result on inter¬ 
changeability of regularized sums and integrals in [LEVE13]. However in contrast 
to the setting considered in [LeVei 3], the individual summands (cu(n, ki,..., k m ) + 
z 2 ) m here depend on the limiting parameter n € N and do not individually admit 
a polyhomogeneous expansion jointly in the summation parameters (ki,k m ) and 
the resolvent parameter z. Only the full sum in the expression of Tr(A n + z 2 )~ m is 
polyhomogeneous, while in the setup of [LEVE13] polyhomogeneity is rather lost 
after summation. 


1.5. Comparison with a theorem by Chinta, Jorgenson and Karlsson. Chinta, Jor¬ 
genson and Karlsson [CJK10] use a different method to establish Theorem 1.6. In 
fact they consider a slightly more general setting where they allow the individual 
cyclic factors Z/uZ to converge to 8 1 at different rates. 

If for rt := (rii,n m ) £ N m we write Tjj 1 = (Z/uiZ) x ■ • ■ x (Z/n m Z), then in our 
picture the resolvent trace of the corresponding combinatorial Laplacian An admits 
a partial polyhomogeneous expansion as (rt, z) —> 00. We consider here the special 
case of Uj = n for all j = 1,m, with the argument for the general setup going 
along the same lines. Then [CJK10] asserts the following. 

Theorem 1.7 ([CJK10]). Consider the rescaled combinatorial Laplacian A^ with cr(A n ) = 
^■a(A^), which approximates the Laplace Beltrami operator A' on a rescaled torus given 
by m copies of f 2 rc)“ 1 8 1 . Then the logarithmic determinant log det A' n admits a regularized 
limit as n —> 00, which equals the logarithm of the zeta-determinant of the A', i.e. 

log deh A' = LIM log det A^. 

n—>00 

This statement can be easily seen to correspond to Theorem 1.6, which is basically 
an issue of conventions. Since the number of non-zero eigenvalues of A n is (n m — 1 ), 
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and by the identity £(s, A') = [In] 2s C(s, A), we obtain the following relations 


log det A^ = log det A n — (n m — 1 ) log 


n 

47 t 2 ' 


logdet^ A' = logdeh A + 2 £( 0 , A) log 27 t. 

Since on closed compact manifolds C( 0 , A) = — dimker A, which equals (— 1 ) in the 
setting of tori, we conclude that Theorems 1.6 and 1.7 are equivalent. 

2. Interchanging regularized limits and integrals 

The main result of this section is presented in [LEVE13, Lemma 3.3], albeit under 
mildly stronger assumptions on the asymptotics. 


Proposition 2.1. Let f e C c 


:,C),R+ = ( 0 ,00), be homogeneous of order d 6 C 


jointly in both variables and asymptotic expansions of the form (1.1) in each of the variables 
z, n — * 00 individually (with the other variable fixed) 

N—1 Mj M 0 

f(z> 1) = Y. Y a i kZ<Xi lo g k M +Y a ok io g k (z) + °( z0cn io g MN M)> 

j=1 k=0 k=0 

(2.l) 

N1 M j 

f(l,n) = Y_ X^ b i knPi lo g k ( n ^ + ^b 0 k log k (n) + o(n pN ' log M *(n)), 

j=1 k=0 k=0 

for some N,N' G N and (cq), ((3j) C C, such that (Re(ctj)), (Re((3j)) are monotonously 
decreasing sequences with Re(a N ) < — l,Re(|3 N /) < min{ 0 , d+ 1 }. Then 


(2.2) 


LIM - f(z, n)dz=- LIMf(z, u)dz +Corr, 


where Corr = J^° f (z, 1 )dz if d = — 1 and zero otherwise. 

Proof Using (2.1) and Re(( 3 N /) < 0 we find for f(z, n) = z d f(l, n/z) 


M' 


LIMf(z, n) = Y (—1 ) k b 0k z d log k (z). 

n—inn » 


k=0 


For any oc G C and k > 0 we compute iteratively 


(2.3) 


L 

z“log k (z)dz = { i=° 

log k 


(—1 ) 3 k! log' J (z) 
(k-j)! (<x + l)i +1 


<x+l 


l Oc+1) ’ 


if a 7^ —1, 


if a = —1. 


Consequently we find 


M' 


( 24 ) 


- LIMf(z, n)dz 


k=0 

0, if d = —1. 
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For the computation of the left hand side in (2.2), we employ the coordinate change 
rule for regularized integrals, cf. [LES97, Lemma 2.1.4], aR d obtain 

'OO poo 

(2.5) ■- f(z,n)dz = n d+1 -- f(x, 1 )dx — a.j 0 ou d+1 log(n) 

Jl Jl/rx 

'OO pi/u 

(2.6) = n d+1 - f(x, 1 )dx — n d+1 -- f(x, 1 )dx — a.j 0 on d+1 log(n), 

.0 Jo 

where <Xj 0 = — 1 . Using (2.1) and (2.3) we arrive for d ^ — 1 at the following expan- 

N'—1 M ( 


sion. 


■1/n 

- f(x, lldx-n^oo 
0 


Y_ X b i knPi_(d+1) Z 

)=1,|3j7d+1 k=0 
N'—1 M j 

Y_ y_ v 

j= 1 ,| 3 j=d +1 k =0 


p, Cd+i) V k! log k ] (V R ) 
j=0 (k — j)! (d + 1 — (3j)’ +1 

.k+l 


log k+1 (1/n) 

(k+l) 


M' 


X b ok n ^Tn 7 — 


k! log M (1/n) 


k=0 


j=0 


(k-j)! (d+1)J +1 


+ o(u | 3n ' (d+1 -log MN (n)). 

If d = — 1 the asymptotic expansion changes slightly to 


■L 


1/ n N1 M j' k 

- f(x, IJdx-n^oo Y_ Y V RPi ~ td+1) 

j=1,|3j4d+1 k=0 

N 1 M j , k+1 M/ , 

+ L Iv ( / ’ 

j=l,(3j=d+l k=0 


k! log k j (1/n) 

^(k-jJKd+i-Pj)^ 1 


M' 


(k+l) 

Jc +1 


y, 1 l k b lo § k+1(1/n) 

k=0 v ’ 

+ o(n |3N ' _(d+1) log M ^ (n)). 

We may now take regularized limit of (2.5) as n —> 00 and find using Re(| 3 N /) < d +1 

M' 

Lt-Lbodd+l U-', ifd#-l, 

if d = — 1. 

Statement follows from comparison of (2.7) and (2.4). □ 


(2.7) 


LIM - f(z, n)dz=< 


k=0 

'OO 

-- f(x, l)dx, 
. 0 
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3 . POLYHOMOGENEOUS EXPANSION OF THE COMBINATORIAL RESOLVENT TRACE 

This section is devoted to a proof of Theorem 1.5, using the Euler Maclaurin 
summation formula in m parameters. Write for u 6 N 0 and x = (xi, ...,x m ) G M m 


cu(n, x) 



Consider any tuple J = ^ of pairwise distinct integers with 

k < m. We write |J| := k and put 

{xj = 0} := {x G N™ fl [0, u] m | Xj, = ... = Xj k = 0}. 

With respect to this notation we obtain for oc G N and z G M 1 the following repre¬ 
sentation of the resolvent trace 


(3.1) Tr(A„+z 2 r" = 2j-l) k 2I ^(cv(n,x) + z 2 r m - 

k=0 |J|=k{x,=0} 

The summand corresponding to k = 0 is given explicitly by 


n rt 

S(z,n) = Y_ Y- >*m) +z 2 ) 

xi =0 x m =0 


We employ the Euler Maclaurin formula for summation in m parameters 2 to derive 
a polyhomogeneous expansion of S(z, n), the other terms in the resolvent trace 
formula (3.1) are treated ad verbatim. For any MgN we obtain 

( 3 - 2 ) S(z,n)= p Pi,i 0 ' ’ • 0 p Pm,m (ai(n,x 1 ,...,x m ) +z 2 ) m , 

(3e{1,2,3,4} m 


where each Pp j acts in the Xj-variable on u G C°°[ 0 > 00) by 


P M U : = 


u(xj)dxj, 


k=W 


(2k)! 

1 


( 2 M + 1 )! 


B 


2M+1 (Xj - [x j ]) 9 ^ 2M+1) u(xj)dx j , 


y (u(xj = n) + u(xj =0)), 

V Z 


if ( 3 j = 1, 

if p, = 2, 

if |3j = 3, 
if (3, = 4. 


Here, B^x) denotes the i-th Bernoulli polynomial and B| the i-th Bernoulli number. 
We study the (z, n)-behaviour of the various summands. Note that for finite u G PI 
all individual operators in the composition Pp^i o ■ ■ • o P|3 miin commute. 


obtained by iterating the standard Euler Maclaurin formula in one single summation parameter. 
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Lemma 3.1. Put f (n, x, z) := tu(n, x) + z 2 . For each keNw may write 

k—1 

1=0 


where H 1)0 ( 3 x .f) 


—<xd x .f and higher functionals H k; e are defined recursively by 
Hk,«( 3 Xj f) = 3 xjHk-i,f ( 3 Xj f) + H k _i £_i ( 3 Xj f) • 3 Xj f, 


where we set H k) k, H k -i = 0 . For k odd we have the estimate 


( 3 - 3 ) 


Hk,f( 3 Xj f)| < Ckf 


n -k+2({+D ? if k >2(£+l), 
x 2(<+ i ) "k j if k<2(£ + 1). 


Proof The recursive structure of 3 xyf _a follows by induction and the only intricate 
statement is the estimate of H k) f (g), g = 3 x .f. Here we introduce a notion of homo¬ 
geneity order for an expression Hk,*(g) by counting for every individual summand in 
H k , £ (g) each additional 3 X . differentiation as lowering homogeneity order by (— 1 ), 
and each factor of g as increasing homogeneity order by (+1). 

With this system, Hi )0 (g) = — ag is of homogeneity order 1 , and by induction 
H k ,{(g) is of homogeneity order (2(£ +1) — k), where in total we count [l + 1) factors 
of g, and (k — (£ + 1 )) derivatives. Note the series expansion 


9 = 



3 Xj g = ^a t (2i 

i=0 


/ 7 tXi \ 2i 


1 ) ( — 

V n 


Consequently, each g in Hk,«(g) carries a factor of Xj, a single additional 3 Xj differ¬ 
entiation applied to g, annihilates that Xj factor, and each further derivative adds a 
factor of ii 1 . Consequently, homogeneity order simply counts the powers of Xj and 
n in the expression for H k ,{(g), leading to the statement. □ 


First, we derive a polyhomogeneous expansion for the summand S(z, n). 
Proposition 3.2. The function S(n, z) admits a partial polyhomogeneous expansion 

m 

S(n,z) = Y h , _ m _ ] {z,n) + H'[z,n), 
j=o 

where each h/e C°° (M 2 ) is homogeneous of order (—m—j) jointly in (z, n), h/(z, 1 ) 
and h/_ m = (1, n) admit an asymptotic expansion of the form (1.1) as z, n —> 00, respectively. 
The remainder term satisfies H^(z,n) = 0 (z^ 2m ^ 2 ), as z —> 00, uniformly in n > 0 . 
Moreover, h/ 2m (z, n) = z^ 2m . 

Proof Since 3 Xj differentiation of odd order applied to (cu(n, x)+z 2 )^ m , always leads 
to factors sin(7rx/u), summands in (3.2) with f 3 j = 2 for some j = 1 ,m vanish. We 
are left to consider summands in (3.2) with | 3 j £ { 1 , 3 , 4 }. Consider first summands 
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in (3.2) with ( 3 j = 3 for some j = 1 ,m. We will use the following basic estimates 


( 3 - 4 ) 

( 3 - 5 ) 


[tu(n, x) + z 2 ) ~xf + z 2 


<c^-L-<c( 2 Z r\ 


[O.n]” 


(tu(n, x)+z 2 ) m d m x < 


n 2 . 2 

/7tXj 

—r sin | 

- 1 

7T 2 

\ n 

2n 

\ m 


-1 


dxi 


y/z 2 + n 2 / 


7T Z 


< ( 27 T) m , 


where we denote all universal constants in the estimates by the same symbol C > 0 . 
Let k = 2 M + 1 > 3 m + 2 be some sufficiently large integer. By Lemma 3.1 (assume 
a > m) we obtain in case k < 2(1 + 1) 


[0,n] 1 


|H v (9 Xj f)f 


-a—1—1 


< 


[0,n] 

C Cz- 2(£+1)+k 


Cx 2 (e+i)-k 

— - (cu(n, x) + z 2 ) 

(cu(n,x ) + z 2 ) Z[l+]) v ' 


2^ —(a+lc—(I+l)) 


[0,n] T1 


(cu(n,x) + z 2 ) 


2 ^-(a+k-(I+1)) 


< Cz k (cu(n, x)+z 2 ) a d m x 
J[0,n] m 


< (3^—E—2(a—m) 


[0,n] r 


(cu(n, x) +z 2 ) m d m x < C( 27 r) m z 


m _—2 a —2 




where we used (3.4) in the second inequality, and employed (3.5) together with 
k > 2 m + 2 in the last inequality In the case k > 2 (£ + 1 ), we compute similarly 


[0,n]" 


|H M (0 x ,f)r*- t " 1 | < Cn- k+2(£+1} (cu(n,x) + z 2 ) 


^ a (£+1) d m x 


2a—2 


< c< 


[0,n ] v 


[0,n]' 


(tu(n, x) + z 2 ) m d m x, if i > m, 


n -k+2(£+1) z -2a-2 


[0,n]’ 


1 d m x, if i < m, 


< Cz 


1 —la.—2 


where we used k > 3 m and (3.5) in the final estimate. Consequently, if | 3 j = 3 for 
some single j = 1 ,m, we may estimate for M > ( 3 m + 1 )/2 uniformly in n > 0 


Ppi.i 0 • ■ ■ 0 P|3 m ,m (oj(n, Xi,x m ) +z 2 ) 


< Cz 


-2a—2 


If #{(3j = 3 } > 1 , the estimates proceed along the same lines. Consider next the 
case with | 3 j = 1 for all j = 1 , ...,m. The corresponding summand is given by the 
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following integral, where we substitute rjj = Xj/u 


Pi.i 


o ■ ■ • o 


l.m 


(o>( 


n, xi, 


x) + z 2 ) 


n 


-2a+m 


[ 0 , 1 ] 


(cu(l,-y) + (z/u) 2 ) 


d m y. 


the expression is homogeneous of order (— 2 a + m) jointly in (n, z). Existence of an 
asymptotic expansion for the homogeneous terms hi 2a+m (-, 1), h/ 2a , m f 1, ■) follows 
e.g. from a careful application of Melrose's push forward theorem 3 . 


It remains to discuss terms with |3j = 4 for some j G { 1 , Each P 4j j is an 

evaluation operator and in fact for any family of pairwise distinct integers J = 
(ji,-,jk) C {1,m} we have 


P 4) j, o--.oP 4)jk (ai(n,x 1 ,...,x m )+z 2 ) “ = (tuj(n,y) + z 2 ) * 
where we have introduced y = (yi, ...,y m _ k ) := (xj k+1 , and write 


cujln 




m—k 


sin" 


i=1 


Ttyi 

n 


Consequently, analysis of the terms with ( 3 j = 4 for j G J proceeds along the lines 
above, with m simply replaced by (m — k). This leads to homogeneous terms of 
homogeneity order (— 2 a + m — k). Setting a = m proves the statement, once we 
observe that the homogeneous term of order (—2m) is given explicitly by 

P 4 ,i o-.-oP 4m (ie(n,x 1 ,...,x m )+z 2 ) m = z~ 2m . 


□ 


We may now prove our first main result. 

Theorem 3.3. The resolvent trace admits a partial polyhomogeneous expansion 

m 

Tr(A n + z 2 r m = Y_ h- m -j(z, n) + H(z, n), 

i=o 

where each h_ m _j G C°° (R 2 ) is homogeneous of order (—m—j) jointly in (z, n), lw m _j (z, 1 ) 
and H_ m _j (1, n) admit an asymptotic expansion of the form (1.1) as z,n-) 00, respectively. 
The remainder term satisfies H N (z, n) = 0 (z -2m ^ 2 ), as z —> 00, uniformly in n > 0 . 
Moreover, h_ 2m (z, n) = 0 . 

Proof. We need to extend the computations in Proposition 3.2 to general terms in the 
expression (3.1). Consider any tuple J = (ji, ...,jk) C { 1 , of pairwise distinct 

integers with k < m and the corresponding term 

Y_ (tu(n,x)+z 2 r m . 

{xj= 0 } 

3 Note that (cufl,y) + (z/nj 2 ) m lifts to a polyhomogeneous function on the blowup space 
[[0,1] m x M + ,A], blown up at the corners A = {(y,t) | y, € (0,1}, t G {0, 00 }}. Pushforward the¬ 
orem of Melrose [MEL 92 , MEL 93 ] then yields an asymptotic expansion of the integral as t —1 0 or 
t —> 00 . The explicit structure of the expansion is irrelevant in our discussion. 
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Put y = (yi,y m _ k ) := (x jk+1 ,x jm ) and note 


Y _ (tu(n,x) +z 2 )~ m = Y . (a>j(n,y) +z 2 )- (m - k) - k 

{x j =0} ye[0,n] m - k 

where we have introduced 


cej(n,y) := 


TXr 


2 m—k 

L» 


sin 


i=1 


Ttyi 


n 


The computations follow for these terms along the lines of Proposition 3.2 with m 
replaced by (m — k) and ex = m. It remains to discuss the homogeneous term of 
order (— 2 m). Note by the statement of Proposition 3.2 on the (— 2 m)-homogeneity 
terms 


-2m(23 n) = Y (-V 


k=0 


L z 

IJMc 


—2m 


L< 

k=0 


m 

k 


2m 


= 0 . 


This proves the statement. 


□ 


4. Convergence of the combinatorial resolvent traces 

In this section we prove Proposition 1.4, using an argument of Dodziuk [Dody6] 
on convergence of zeta functions. 

Proposition 4.1. For any integer <x > m 

lim 0 < k) Tr(A n + z 2 )~ a = 0 < k) Tr(A + z 2 )“ a . 

TL—>00 

Proof. The argument does not depend on the explicit structure of the spectra for A n 
and A. Hence we write 


o-(A n )={A k , n |k = 0 ,...,N(n)}. 
a(A)={A k |keN 0 }, 

where both sets are ordered in an ascending order, N (n) is the total number of A n - 
eigenvalues counted with their multiplicities, N (n) —> 00 monotonously increasing 
and A kj n —> A k as n —> 00, for each fixed k e No. 

Fix any e > 0 . Since the resolvent trace 

OO 

Tr(A + z 2 )- K = ^(A k + z 2 )- a , 

k=0 

is a convergent series, there exists K e N sufficiently large, such that 

OO 

Y_[\ + z 2 ) “ < e/3. 

k=K 
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Consequently, since (A kin ) neN is an ascending sequence, by the minimax principle 
(cf. [Dod76]), we may estimate 

OO N (n) 

0<^(A k + z 2 r a - X.( Ak,n + z 2 r a < 2 e/ 3 . 

k=K k=K 

Finally, given convergence of the spectrum, choose n > n 0 , such that 

K 

0<^((A k + z 2 r a - (\n + z 2 r a ) < e/3. 

k=0 

Thus, for any given e > 0 , there exists u £ N sufficiently large, such that 

00 N( n ) 

0 < Tr(A + z 2 r a - Tr(A n + z 2 r a = ^ 1 (A k + z 2 )~ a - ^ (V + ^ < e. 

k=0 k=0 

Hence the combinatorial resolvent trace Tr(A n + z 2 )~ a converges to Tr (A + z 2 ) ~ a as 
u —> 00. The general statement follows from the observation 

3 z Tr(A n + z 2 )“ a = 2 zTr(A n + z 2 )-^ 1 , 

3 z Tr(A + z 2 )~ a = 2 zTr(A + z 2 )^ 1 • 

□ 


5. Proof of the main result 
We may now prove our main theorem. 

Theorem 5.1. The logarithmic determinant logdet A n admits a regularized limit as n 
00 and 

log deh A = LIM log det A n . 

n —>00 

Proof. By (1.6) and Proposition 1.4 


log det c A = —2 j- z 2m-1 Tr(A + z 2 )- m dz 


,2^—m. 


0 

poo 


— — 2 — z 


,2m—1 


lim Tr(A n + z 2 )' m dz 


0 

'OO 


= -2 - z 2m ^‘ LIMTr(A n + z 2 )“ m dz. 

q n—»oo 

By Theorem 1.5, Tr(A n + z 2 )' m admits a polyhomogeneous expansion 

m 

Tr(A n + z 2 r m = Y h_ m _j(z, n) + H(z, n), 

j=0 
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where each h m j G C°°(R 2 ) is homogeneous of order (—m — j) jointly in (z, n), and 
the remainder term satisfies H(z, n) = 0 (z 2m_2 ) as z —> oo, uniformly in n > 0 . By 
Proposition 1.3 


-2 - z 2 - 1 LIM V h_ m _j (z, n) dz = — 2 LIM 

, n—inn * n— 


j=o 


00 ™ 

2 m— 


h z' 

1 


1 Y h_ m _j(z,n)dz 
j=o 


-2-- z 


,2m—1 


h—2m (z, 1 ) dz. 


Each homogeneous term h m .j admits a regularized limit as n —> 00 and hence, 
by Proposition 1.4, the remainder term H(z, n) admits a regularized limit as n —> 00 
as well. Since the estimate H(z, n) = 0 (z~' 2m_2 ) for z —> 00 is uniform in n > 0 , the 
regularized limit LIM n ^oo H(z, n) is in fact a true limit lim n _ )00 H(z, n). We find by 
dominated convergence 


-2 


z 2m 1 lim H(z,n)dz = —2 lim 

TL—>00 Tl—>00 


z 2m ^(z, n)dz. 


Consequenty, we arrive at the following intermediate result 


- 2 f z 2m_1 Tr(A + z 2 )~ m dz = — 2 - 1 - z 2 " 1 - 1 ■ ~ 2 ^ m 


,2m— 1 


lim Tr(A n + z j m dz 

TL—>00 


( 5 - 1 ) 


= - 2 LIM f z 2m_1 Tr(A n + z 2 )^ m dz 


2\—m, 


- 1 -- Z' 


,2m—11 


-2m 


Z, l)dz. 


Note that TrfA + z 2 ) m ~ z 2m as z —> 0 and hence z 2m lr TrfA + z 2 )‘ m is not integrable 
at z = 0 due to non-trivial kernel of A. Similar statement holds for Tr(A n + z 2 )' m . 
Subtracting the contribution from harmonic functions 4 , we find 


f 1 

- 2 - z' 

Jo 

rl 


2m—1 


Tr(A + z 2 )^ m dz = -2 


,2m—1 


(Tr(A 


z 2 r m -z 


-2m 


) dz, 


—2 - z 2m-1 Tr(A n + z 2 )^ m dz = -2 

Jo 


Z 2 m—1 ( Tr(An + z 2)-m _ z -2m) dz? 


where the integrals on the right hand side exist in the usual sense. Consequently, re¬ 
placing regularized integrals with the usual integrals, we may interchange integrals 
and limits. This yields 


( 5 - 2 ) 


—2 - z Zm_ Tr(A + zT m dz = - 2 - z /m ~' lim Tr(A n + z 


2 3 —m 


dz 


= - 2 lim f- z"Tr(A n + z 2 )~ m dz. 

n—>00 


4 Note from (1.3) that ker A n and ker A are both one-dimensional. 
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By (5.1) and (5.2), we find in view of the formula (1.6) 


( 5 - 3 ) 


log det^ A = LIM log det A n 

C Tl—>00 



Z 2m l)dz. 


The statement now follows from the fact that h 2m(z, n) = 0 by Theorem 3.3. 


□ 
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